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Near-infrared  diode  laser  absorption  diagnostic 
for  temperature  and  water  vapor  in  a 
scramjet  combustor 


Jonathan  T.  C.  Liu,  Gregory  B.  Rieker,  Jay  B.  Jeffries,  Mark  R.  Gruber, 
Campbell  D.  Carter,  Tarun  Mathur,  and  Ronald  K.  Hanson 


Tunable  diode  laser  absorption  measurements  of  gas  temperature  and  water  concentration  were  made  at 
the  exit  of  a  model  scramjet  combustor  fueled  on  JP-7.  Multiplexed,  fiber-coupled,  near-infrared  distrib¬ 
uted  feedback  lasers  were  used  to  probe  three  water  vapor  absorption  features  in  the  1.34-1.47  |jum 
spectral  region  (2v1  and  vx  +  v3  overtone  bands).  Ratio  thermometry  was  performed  using  direct- 
absorption  wavelength  scans  of  isolated  features  at  a  4-kHz  repetition  rate,  as  well  as  2 f  wavelength 
modulation  scans  at  a  2-kHz  scan  rate.  Large  signal-to-noise  ratios  demonstrate  the  ability  of  the 
optimally  engineered  optical  hardware  to  reject  beam  steering  and  vibration  noise.  Successful  measure¬ 
ments  were  made  at  full  combustion  conditions  for  a  variety  of  fuel/air  equivalence  ratios  and  at  eight 
vertical  positions  in  the  duct  to  investigate  spatial  uniformity.  The  use  of  three  water  vapor  absorption 
features  allowed  for  preliminary  estimates  of  temperature  distributions  along  the  line  of  sight.  The 
improved  signal  quality  afforded  by  2 f  measurements,  in  the  case  of  weak  absorption,  demonstrates  the 
utility  of  a  scanned  wavelength  modulation  strategy  in  such  situations.  ©  2005  Optical  Society  of 
America 

OCIS  codes:  120.0120,  120.1740,  300.0300,  300.1030,  300.6260,  300.6380. 


1.  Introduction 

The  history  of  scramjet  engine  research  spans  more 
than  half  a  century  and  continues  to  generate  inter¬ 
est  as  an  enabling  technology  for  air-breathing  hy¬ 
personic  propulsion.  Although  the  basic  concept  of  a 
scramjet  engine  is  straightforward,  research  in  this 
field  has  been  hampered  by  a  number  of  challenges. 
Coupled  interactions  between  shocks  and  boundary 
layers  exist  within  the  engine,  such  that  complicated 
engine  geometries  are  required  for  robust  flame  hold¬ 
ing  and  heat  release.  A  recent  review  of  scramjet 
engine  development  outlines  some  of  these  issues.1 

Because  of  the  many  challenges  and  unknowns  in 
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scramjet  engine  development,  diagnostic  tools  are  in¬ 
valuable  for  expediting  the  design  and  testing  pro¬ 
cess.  For  any  combustion  device,  gas  temperature  is 
paramount  for  diagnosing  combustion  efficiency.  Wa¬ 
ter  vapor  is  a  major  combustion  product  that  is  also  a 
key  indicator  of  the  extent  of  combustion.  This  paper 
describes  a  diode  laser  absorption  sensor  that  inter¬ 
rogates  multiple  water  vapor  transitions  in  the  near- 
infrared  spectral  region.  From  the  ratio  of  the 
measured  line  strengths  of  individual  spectral  fea¬ 
tures,  path-integrated  temperatures  may  be  inferred, 
along  with  water  vapor  concentrations  in  the  post¬ 
combustion  region  of  a  scramjet  test  facility  located  at 
Wright-Patterson  Air  Force  Base  (WPAFB)  in  Ohio. 

Diode  laser  sensors  have  been  gaining  popularity 
in  a  wide  array  of  practical  measurement  applica¬ 
tions  ranging  from  trace-gas  detection2-4  and  com¬ 
bustion  diagnostics5-13  to  numerous  biomedical 
applications.  Their  low  cost,  small  size,  and  simple 
turn-key  operation  enable  the  development  of  optical 
sensors  that  are  feasible  for  a  wide  range  of  potential 
applications.  Fiber-coupled  diode  lasers  are  espe¬ 
cially  convenient  to  route  laser  radiation  in  and 
around  test  facilities.  For  probing  water  vapor  fea¬ 
tures,  telecommunications  distributed  feedback 
(DFB)  lasers  exist  in  the  1.3  to  1.5-|xm  region,  pro- 
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Fig.  1.  Left,  detected  signal  for  a  direct-absorption  scan  (one  scan)  near  1392  nm.  (b)  Right,  a  typical  absorbance  plot  for  an  isolated 
absorption  feature. 


viding  access  to  water  vapor  overtone  and  combina¬ 
tion  bands  (2iq  and  u1  +  v3).  These  DFB  lasers  are 
narrowband,  single  mode,  and  relatively  powerful 
(10-30  mW)  and  may  be  tuned  over  several  wave- 
numbers  at  kilohertz  rates  by  modulating  the  injec¬ 
tion  current.  Static  pressures  in  the  postcombustion 
measurement  location  in  a  scramjet  engine  are  on  the 
order  of  1  atm,  where  the  linewidths  of  water  absorp¬ 
tion  features  are  of  order  0.1  cm-1.  Therefore,  wave¬ 
length  scans  across  absorption  features  at  kilohertz 
rates  can  be  easily  accomplished  with  the  DFB  lasers 
used  in  this  study. 

Many  challenges  complicate  the  implementation  of 
laser-based  sensors  in  realistic  engine  facilities  and 
are  dealt  with  in  this  work: 

1.  The  turbulent  flow  fields  of  many  combustors, 
such  as  the  scramjet,  in  which  large  boundary  layers 
and  shock  structures  exist,  lead  to  variable  density 
gradients  within  the  flow  that  cause  beam  steering  of 
probe  lasers.  Mechanical  vibrations  may  cause  addi¬ 
tional  deviations  in  laser  beam  paths.  However,  the 
optical  hardware  can  be  designed  to  reject  such  noise 
sources  and  enable  the  acquisition  of  data  of  high 
signal-to-noise  ratio  (SNR),  as  described  in  Subsec¬ 
tion  3. A  and  demonstrated  in  Subsection  4. A. 

2.  Limited  optical  access  and  difficult  mounting 
environments  make  fiber-optic  transport  of  laser  ra¬ 
diation  extremely  valuable.  Single-mode  fibers  are 
favored  for  laser  beam  transmission,  as  a  spatially 
smooth  Gaussian  output  is  desired.  However,  multi- 
mode  fibers  are  favored  for  light  collection  in  order  to 
provide  a  larger  target  for  the  collection  of  light  in  the 
presence  of  various  beam-steering  and  misalignment 
mechanisms.  Unfortunately,  the  light  beam  exiting  a 
multimode  fiber  exhibits  mode  noise,  but  this  can  be 
effectively  suppressed  by  eliminating  any  losses  after 
the  transmitted  probe  beam  is  collected  (see  Subsec¬ 
tion  3.  A). 

3.  Temperature  and  concentration  gradients  due, 
for  example,  to  boundary  layers,  can  complicate  the 
interpretation  of  line-of-sight  absorption  ratios  for 
gas  temperatures.  The  use  of  multiple  lines  provides 
additional  information  about  the  uniformity  of  tem- 
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perature  distributions  along  the  line  of  sight,  as  will 
be  discussed  in  Subsection  4.C. 

4.  The  SNR  in  scanned  direct  absorption  mea¬ 
surements  is  often  limited  by  errors  associated  with 
the  baseline  fitting  of  absorption  features.  Derivative 
techniques  like  2 f  spectroscopy  offer  benefits  in  cases 
of  weak  absorption,14  as  described  in  Subsection  4.D. 

The  proof  of  principle  of  this  sensor  is  shown  in 
Subsection  4.B,  where  example  results  are  shown  of 
temperature  and  water  concentration  time  histories 
for  various  combustor  conditions  and  combustor  lo¬ 
cations.  The  good  temporal  resolution  of  these  mea¬ 
surements  allows  for  the  frequency  power  spectrum 
analysis  of  the  stability  of  the  gas  flow. 

2.  Measurement  Principles 

A.  Direct  Absorption  Methods 

Scanned-wavelength,  direct-absorption  sensors  are 
useful  in  measurement  applications  in  which  isolated 
and  narrow  spectral  features  exist.  A  linear  variation 
in  injection  current  is  used  to  scan  DFB  diode  lasers 
simultaneously  in  intensity  and  wavelength.  Figure 
1(a)  shows  a  typical  direct-absorption  scan  using  a 
diode  laser,  in  which  the  laser  intensity  changes  in 
response  to  an  injection  current  ramp.  A  polynomial 
fit  to  the  nonabsorbing  wings  of  the  absorption  fea¬ 
ture  is  used  to  extrapolate  a  zero-absorption  baseline, 
from  which  an  absorbance  plot  of  the  absorption  fea¬ 
ture  may  be  constructed,  as  shown  in  Fig.  1. 

The  fundamental  equation  for  absorption  spectros¬ 
copy  is  the  Beer-Lambert  Law.  The  Beer-Lambert 
relation  gives  the  transmitted  intensity,  It ,  of  mono¬ 
chromatic  radiation  after  passing  through  an  absorb¬ 
ing  gas,  given  here  for  a  uniform  gas  medium, 

T(i>)  =  (/”)  =  exp[-a(u)L],  (1) 

where  t(v)  is  the  transmission  coefficient,  I0  is  the 
incident  radiation  intensity,  a(u)  is  the  absorption 
coefficient,  L  is  the  optical  absorption  path  length, 
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and  absorbance  is  defined  as  the  quantity  a (v)L.  The 
absorption  coefficient  is  the  product  of  line  strength, 
S[cm-2/atm],  absorbing  species  partial  pressure, 
P;[atm],  and  line-shape  function  4>(u)  [cm]:  a(v) 
=  S$(y)pi.  The  absorption  coefficient  is  a  function  of 
temperature,  pressure,  and  gas  composition,  which 
enables  these  parameters  to  be  inferred  through  care¬ 
fully  designed  spectroscopic  measurements. 

The  line  strength  of  an  absorption  feature  is  a  func¬ 
tion  of  temperature;  the  transition  frequency,  v0;  the 
lower-state  energy  of  the  quantum  transition,  E "  as 
well  as  the  partition  function  of  the  species,  Q(T) 
(rovibrational  states  are  relevant  for  near-infrared 
water  vapor  spectra): 


S(T)  S(T0 )  T 


hcv  o 

T0  Q(T0)  expl  KT 


1  -  exp  - 


hcv0 

KTo 


X  exp 


he 

'7f e"[t~Yo 


(2) 


Here,  K  is  Boltzmann’s  constant,  S(T0)  is  the  line 
strength  at  a  reference  temperature  T0  (generally 
tabulated  at  room  temperature,  296  K),  h  is  Planck’s 
constant,  and  c  is  the  speed  of  light.  The  partition 
function  has  been  analytically  determined  by  various 
researchers  who  have  provided  correlations  used  in 
our  spectral  simulations.15 16  For  water  vapor  transi¬ 
tions,  the  HITRAN  database17’18  provides  a  compila¬ 
tion  of  spectroscopic  parameters  such  as  line  strength 
and  broadening  coefficients.  However,  in  our  experi¬ 
ence,  an  accurate  temperature  measurement  re¬ 
quires  the  validation  of  the  HITRAN  parameters  for 
selected  transitions.19 

In  practical  measurements,  the  spectrally  inte¬ 
grated  absorbance  (area)  of  a  transition  is  often  mea¬ 
sured,  which  is  directly  proportional  to  line  strength: 


area  = 


|di>  S  (T)pSpecjesZ/. 


(3) 


The  ratio  of  the  integrated  absorbance  of  two  transi¬ 
tions  is  only  a  function  of  temperature  for  multi¬ 
plexed  laser  beams  probing  identical  paths: 

areax  S^T) 

- =  — — —  (4) 

area2  S^T) 

By  measuring  the  integrated  absorbance  of  two  tran¬ 
sitions,  temperature  is  inferred  with  Eq.  (4),  and  spe¬ 
cies  concentration  is  inferred  with  Eq.  (3). 

B.  2 f  Methods 

Wavelength  modulation  spectroscopy  (WMS)  with  2 f 
detection  is  a  common  technique  used  to  improve 
SNR  in  measurements  as  it  shifts  the  detection  band 
to  higher  frequencies,  is  phase  sensitive,  and  pro¬ 


vides  a  zero  baseline  sensing  strategy.  However,  the 
improved  SNR  comes  at  the  expense  of  additional 
requirements  for  calibration.  The  scramjet  engine  is 
well  suited  for  utilizing  scanned-2/*  ratio  thermome¬ 
try  because  static  pressures  are  relatively  low  and 
absorption  features  are  consequently  narrow  and  iso¬ 
lated.  The  theory  of  WMS  will  not  be  described  in 
detail  here,  as  it  is  well  documented  in  the 
literature,20-28  along  with  various  practical  imple¬ 
mentations  of  WMS  for  gas  sensing  using  diode 
lasers.29-35 

Since  the  scramjet  facility  investigated  in  our  re¬ 
search  provided  a  long  measurement  path  length 
(23  cm),  near-atmospheric  pressures,  large  water 
vapor  mole  fractions  (—10%),  and  moderate  temper¬ 
atures  (<2000  K),  absorption  strengths  were  gener¬ 
ally  large  enough  that  direct-absorption  methods 
were  ideal.  Although  direct-absorption  methods  are 
preferred  for  the  two  strongest  lines  used  in  this 
study,  analysis  of  the  weakest  line  demonstrates  the 
utility  of  2 f  methods  for  improving  SNR  for  low- 
absorbance  measurements. 

3.  Experimental 

A.  Optical  Hardware 

For  scanned-wavelength  temperature  and  concentra¬ 
tion  sensors,  fiber-coupled,  DFB  diode  lasers  from 
NEL  America  Inc.  are  multiplexed  using  a  standard 
single-mode  fiber  combiner  (1x4)  designed  for  the 
near  infrared.  Three  of  the  lasers  (1392,  1343, 
and  1469  nm)  are  temperature  tuned  to  lase  at  wave¬ 
lengths  near  selected  absorption  features.  An  exter¬ 
nal  modulation,  consisting  of  either  a  4-kHz  saw 
tooth  ramp,  for  direct-absorption  scans,  or  a  2-kHz 
ramp  summed  with  a  faster  190-kHz  sinusoidal 
modulation,  for  scanned-2/*  spectroscopy,  is  fed  into 
each  of  the  current-source  units.  A  fourth  laser  is  a 
nonresonant  laser  near  1315  nm  that  is  used  to  mon¬ 
itor  transmission  fluctuations. 

Fiber  optics  are  useful  for  routing  laser  radiation 
with  minimal  intrusion  such  that  transmission  and 
collection  optics  may  be  easily  mounted  on  the  engine 
test  stand.  A  nonfibered  approach  using  tree-space 
optics  is  not  feasible  for  most  realistic  engine  diag¬ 
nostic  applications.  The  fiber-coupled  (9-pm  single¬ 
mode  fiber)  multiplexed  beam  is  transmitted  across  a 
test  region  using  an  aspheric  collimator  (Thorlabs 
F230FC-C).  On  the  collection  side,  a  larger  collimator 
(Oz  Optics  HPUCO-25-1300-M-10BQ)  is  used  to  focus 
the  free-space  beam  into  a  400- jxm  diameter  multi- 
mode  fiber  with  a  numerical  aperture  (NA)  of  0.39 
(Thorlabs  M20L05).  The  1-cm  diameter,  1-cm  focal 
length,  Oz  Optics  collection  lens  offers  a  large  “sweet 
spot”  for  the  collection  of  the  collimated  laser  beam, 
such  that  the  collection  is  tolerant  of  relatively  large 
changes  in  lateral  and  angular  alignment.  This  is  a 
key  component  for  minimizing  beam-steering- 
induced  laser  intensity  fluctuations  by  the  turbulent 
flow  field  within  the  combustor.  A  400- jxm  fiber 
routes  the  collected  signal  into  a  grating-based  de¬ 
multiplexing  setup,  where  a  large  aspheric  lens  (Op- 
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Fig.  2.  Optical  configuration  for  two  multiplexed  lasers.  The  ac¬ 
tual  sensor  used  in  these  measurements  contained  four  multi¬ 
plexed  lasers. 


Expanded  wet  air, 


Fig.  3.  Block  diagram  of  the  SCRAM  JET  facility. 


tosigma  023-2392)  is  used  to  collimate  the  output 
from  the  400-pm  fiber  and  send  it  onto  a  30-mm 
square  diffraction  grating  (first  order)  with  1200 
grooves/mm  (Edmund  Optics  NT43-852).  The  de¬ 
multiplexed  beams  are  then  dispersed  in  space 
and  focused  with  concave  mirrors  onto  individual 
4-MHz  large-area  (3-mm  diameter)  amplified  In- 
GaAs  photodetectors  from  Electro-Optical  Systems 
Inc.  (IGA-030-E5/ 4MHz). 

A  benefit  of  a  grating-based  dispersion  system  is 
that  it  allows  for  the  rejection  of  broadband  emission 


from  a  combustor.  The  dispersion  system  behaves 
approximately  as  a  10-nm  bandpass  filter  to  reject 
any  flame  luminosity  captured  by  the  collection  fiber. 
The  entire  dispersion  apparatus  is  enclosed  in  a  tent 
and  purged  with  a  nitrogen  flow  to  remove  ambient 
water  vapor  from  the  dispersion  beam  paths.  Figure 
2  illustrates  the  experimental  configuration. 

The  high-bandwidth  large-area  detectors  insure 
collection  of  the  entire  demultiplexed  beam,  which 
minimizes  fiber  mode  noise  introduced  by  the 
400- fxm  collection  fiber.  Mode  noise  results  from  the 
speckled  structure  of  light  emitted  from  multimode 
fibers  owing  to  interference  between  rays  of  light 
traveling  in  different  cavity  modes  within  the  fiber. 
Since  this  structure  is  wavelength  dependent,  it  can 
lead  to  intensity  distortions  in  diode  laser  wavelength 
scans  if  only  a  portion  of  the  light  emitted  from  a 
multimode  fiber  is  detected. 

Detector  signals  from  all  channels  are  simulta¬ 
neously  recorded  at  a  4-MHz  sampling  rate  by  a  Na¬ 
tional  Instruments  data-acquisition  system  using  a 
Lab  VIEW  scope  program.  The  2/*  line  shapes  are  re¬ 
covered  by  running  a  digital  lock-in  program  on  Lab- 
VIEW  with  a  low-pass  filter  time  constant  of  70  kHz. 

B.  Scramjet  Facility 

A  detailed  description  of  the  scramjet  facility  at 
WPAFB  will  not  be  presented  here.  A  more  com¬ 
plete  account  of  the  facility  may  be  found  in  the 
literature.36"39  A  block  diagram  and  schematic  of 
the  scramjet  facility  at  WPAFB  is  shown  in  Figs.  3 
and  4,  respectively.  An  electrically  heated  dry 
compressed-air  flow  is  continuously  supplied  to  the 
direct-connect  combustion  facility  at  13.6  kg/s,  up 
to  5.2  MPa  and  920  K  total  pressure  and  tempera¬ 
ture.  The  flow  is  further  heated  by  combustion  in 
a  vitiator,  which  results  in  a  high-temperature 
(up  to  2500  K)  and  high-pressure  flow.  Oxygen  is 
added  to  replenish  that  which  was  burned  in  the 
vitiator.  The  hot,  high-pressure  flow  is  then  expanded 
to  supersonic  speeds  where  it  enters  a  model  scramjet 
combustor.  The  JP-7  fuel  is  added  and  burned  in  a 
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H- 3.95”  (10.03  cm) 
G  -  3.70”  (9.40  cm) 

F  -  3.325”  (8.45  cm) 

E  -  2.575"  (6.54  cm) 

D- 1.575”  (4.00  cm) 
C- 0.825"  (2.10  cm) 
B  -  0.45"  (1.14  cm) 
A -0.20”  (0.51  cm) 


Fig.  5.  Location  of  eight  measurement  beam  paths  through  the 
optical  access  windows  (vertical  positions  A  through  H). 


cavity-based  flame  holder.  Engine  runs  for  the  mea¬ 
surements  described  in  this  work  were  of  several  sec¬ 
onds  in  duration. 

The  postcombustion  gases  are  monitored  down¬ 
stream  via  a  rectangular  2.5  cm  X  12  cm  quartz  win¬ 
dow.  Collimation  optics  for  the  transmission  and 
collection  of  laser  radiation  were  offset  from  the  com¬ 
bustor  windows  by  at  least  6  in.  to  avoid  damage  due 
to  radiation  from  flame  luminosity.  Nitrogen  purge 
tubes  were  installed  to  eliminate  atmospheric  water 
interference  in  the  sensor’s  free-space  laser  paths. 
Measurements  were  made  along  one  of  eight  horizon¬ 
tal  spanwise  beam  paths  positioned  vertically  along 
the  rectangular  optical  access  region,  as  shown  in 
Fig.  5.  Transmission  and  collection  optics  were 
mounted  on  vertical  translation  stages  on  either  side 
of  the  combustor,  so  that  beam  paths  could  be  moved 
between  combustor  runs.  The  precision  of  the  trans¬ 
lation  stages  was  sufficient  to  maintain  alignment  at 
all  measurement  locations. 

C.  Line  Selection  and  Validation 

Spectral  features  must  be  selected  possessing  a  com¬ 
bination  of  lower-state  energies  such  that  the  ratio  of 
recorded  signal  strengths  is  a  sensitive  function  of 
temperature.  Accurate  determination  of  absorption 
baselines  requires  isolated  spectral  features.  A  dis¬ 
cussion  of  the  strategy  behind  this  optimal  selection 
procedure  may  be  found  in  the  literature.5’40  41 

The  spectral  features  used  in  this  sensor  are  tab¬ 
ulated  in  Table  1.  Lower-state  energies,  E ",  are  given 
by  HITRAN04.17 18  Line  strengths  are  tabulated  from 
HITRAN04,  as  well  as  from  our  own  laboratory  val¬ 


Temperature  (K) 

Fig.  6.  Line  strength  versus  temperature  for  water  features  at 
1392,  1343,  and  1469  nm,  using  validated  parameters  (see  Table 
1). 


idations.7’11  A  plot  of  the  line  strength  versus  tem¬ 
perature  for  the  three  water  features  used  in  this 
study,  with  laboratory-validated  parameters,  is 
shown  in  Fig.  6. 

4.  Results 

A.  Improvements  to  Signal-to-Noise  Ratio 

The  water  vapor  absorption  data  collected  here  has 
significantly  improved  SNR  compared  with  earlier 
measurements42  because  of  improved  optical  engineer¬ 
ing.  Collection  optics  had  a  large  sweet  spot  for  cap¬ 
turing  the  collimated  laser  beam,  resulting  in  excellent 
tolerance  to  beam  steering.  Figure  7  shows  a  compar¬ 
ison  between  the  single-scan  transmission  signals  re¬ 
ported  here  (left)  and  the  signals  from  a  previous 
measurement  campaign  (right).42  The  current  SNR  is 
approximately  35  times  larger  than  in  previous  data. 
The  scan  rates  of  the  current  measurements  [Fig.  7, 
right]  are  also  four  times  faster  (4  kHz)  than  the  scan 
rates  utilized  in  the  earlier  measurements  [Fig.  7, 
left].  The  most  important  changes  are  the  use  of  a 
larger-diameter  collection  lens  (1  cm),  a  larger- 
diameter  collection  fiber  (400- |xm  versus  50- jxm), 
and  a  collection  fiber  with  a  larger  NA  (0.39  versus 
0.22). 


Table  1.  Line  Strength  and  Lower-State  Energy  Parameters  for  Water  Vapor  Features  Used  in  Measurement  Campaigns 


Wavelength 

(nm) 

HITRAN04 

Wavenumber 

(cm^1) 

HITRAN04 

Line  Strength,® 

S  (cm^2  atm-1) 
HITRAN04 

Line  Strength,® 

S  (cm  2  atm-1) 
Measured 

Lower-State 

Energy  (cm1) 
HITRAN04 

1392 

7185.60 

1.97  (lor2 

1.96  (10)“2 

1045.1 

1343 

7444.35+ 

i.i2  cion3 

1.10  (10)“3 

1774.8 

7444.37  (combined) 

1806.7 

1469 

6807.83 

1.02  (10)-® 

6.48  (lor7 

3319.4 

a296  K. 
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Fig.  7.  Left,  single  multiplexed  absorption  scan  taken  at  1  kHz  (1-ms 
taken  at  4  kHz  (0.25-ms  scan).  Data  taken  September  2004. 


B.  Engine  Diagnostics 

Temperature  and  water  concentration  time  histories, 
at  a  4-kHz  bandwidth,  were  acquired  to  characterize 
the  vitiator,  as  well  as  for  full  combustion  runs  in  the 
scramjet  facility.  An  example  plot  of  measured  tem¬ 
perature  and  water  partial  pressure  in  a  vitiator-only 
flow  is  shown  in  Fig.  8.  The  measurement  beam  path 
was  at  position  C  (see  Fig.  5),  and  the  average  static 
pressure  at  the  window,  measured  with  pressure 
taps,  was  0.22  ±  0.01  atm.  Water  vapor  features  at 
1392  and  1343  nm  were  used  to  infer  temperature, 
as  they  are  relatively  strong  at  these  temperatures 
and  exhibit  good  temperature  sensitivity.  Tempera¬ 
tures  were  also  inferred  by  performing  Voigt  fits  to 
the  absorption  features  and  using  absorbance  area 
(spectrally  integrated  absorbance)  ratios,  which  re¬ 
moves  the  need  to  know  total  pressure  or  gas  compo¬ 
sition.  Absorption  peak-height  ratios  were  also  used 
to  infer  temperatures,  utilizing  the  spectroscopic  pa¬ 
rameters  listed  in  HITRAN04.  Both  the  peak-height 
ratio  and  the  area  ratio  data-analysis  methods  yield 
similar  temperature  results  for  the  vitiator-only  mea¬ 
surements.  Water  partial  pressure  measurements 


Fig.  8.  Vitiator  only:  Temperature  versus  time  (left  axis)  and 
water  partial  pressure  versus  time  (right  axis);  4-kHz  bandwidth; 
position  C. 
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scan).  Data  taken  9/2002.  Right,  single  multiplexed  absorption  scan 


may  be  inferred  from  temperature  measurements 
and  the  integrated  absorbance  of  the  Voigt-fitted  wa¬ 
ter  features. 

For  full-combustion  runs,  Voigt-fitted  area  ratios  are 
used  to  infer  temperature.  An  example  of  a  tempera¬ 
ture  and  water  partial  pressure  time  history  measure¬ 
ment  is  shown  in  Fig.  9.  These  measurements  utilize 
the  absorption  features  at  1392  and  1343  nm.  A 
1-kHz  filter  is  used  to  smooth  the  data  recorded  at  a 
4-kHz  bandwidth. 

Temperature  and  water  partial  pressure  measure¬ 
ments  were  made  across  a  fixed  beam  path  at  position 
C  (see  Fig.  5),  as  a  function  of  combustor  fuel/air 
equivalence  ratio.  Average  results  are  plotted  in  Fig. 
10,  demonstrating  the  ability  of  the  sensor  to  capture 
varying  engine  conditions,  through  changes  in  both 
temperature  and  water  concentration.  Expected  tem¬ 
peratures  and  water  concentrations,  based  on  a  ID 
combustion  model,  are  shown  for  comparison.  These 
simulations  are  calculated  by  scramjet  researchers  at 
WPAFB,  utilizing  information  from  various  engine 
sensors,  such  as  pressure  taps  and  flow  rate  moni- 


0-| - , - j - 1 - ( - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 0,00 

0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5 


Time  (s) 

Fig.  9.  Full  combustion  run:  Temperature  versus  time  (left  axis) 
and  water  partial  pressure  versus  time  (right  axis).  Fuel/air  equiv¬ 
alence  ratio:  =  1.02;  position  C;  4-kHz  measurement  bandwidth 
filtered  to  1  kHz  in  this  plot. 
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Fig.  10.  Left,  average  temperature  versus  combustor  fuel/air  equivalence  ratio,  <F.  Measured  and  expected  results  (based  on  a  combustor 
ID  model)  are  shown.  The  error  bars  represent  a  standard  deviation  in  the  data  at  the  full  4-kHz  bandwidth;  position  C.  Right,  average 
water  partial  pressure  versus  combustor  fuel/air  equivalence  ratio,  <F.  Measured  and  expected  results  (based  on  a  combustor  ID  model) 
are  shown.  The  error  bars  represent  a  standard  deviation  in  the  data  at  the  full  4-kHz  bandwidth;  position  C. 


tors.  There  is  a  slight  discrepancy  between  measured 
and  expected  results  at  the  lower  equivalence  ratios, 
which  should  be  a  subject  for  further  investigation. 
Temperature  and  water  concentrations  were  also 
measured  across  the  scramjet  combustor  at  the  eight 
vertical  locations  shown  in  Fig.  5.  Examples  are 
shown  in  Fig.  11  with  error  bars  to  indicate  measure¬ 
ment  uncertainty  (one  standard  deviation  at  the  full 
4-kHz  bandwidth).  The  presence  of  larger  density 
gradients  at  the  top  of  the  combustor  causes  signal 
quality  to  degrade  because  of  beam-steering  noise, 
especially  at  the  upper  two  measurement  locations. 

The  4-kHz  measurement  bandwidth  of  the  direct- 
absorption  temperature  measurements  allows  for  an 
analysis  of  frequency  content  of  up  to  2  kHz.  A  fast 
fourier  transform  (FFT)  power  spectrum  of  the  data 
shown  in  Fig.  9  indicates  low-frequency  modes  below 
150  Hz,  suggesting  mechanical  vibration  frequencies 
or  combustion  instabilities.  An  example  power  spec¬ 
trum  is  shown  in  Fig.  12,  where  distinct  frequency 
modes  are  noticeable  at  —21,  —63,  and —116  Hz. 


Vertical  position  from  bottom  (ins) 

Fig.  11.  Time-averaged  temperature  (left  axis)  and  water  partial 
pressure  (right  axis)  versus  vertical  measurement  path  location. 
Fuel/air  equivalence  ratio:  ~  1.00;  4-kHz;  bandwidth. 


Such  frequency  content  is  of  value  to  engine  design¬ 
ers  who  desire  to  suppress  such  fluctuations  by  me¬ 
chanical  design  and/or  active  combustion  control. 
These  results  illustrate  the  value  of  time-resolved 
temperature  measurements  for  aerospace  ground 
test  experiments. 

C.  Line-of-Sight  Temperature  Distribution 

The  scramjet  combustor,  as  with  nearly  all  practical 
combustors,  does  not  exhibit  a  uniform  temperature 
profile  across  the  measurement  path.  The  vertically 
scanned  temperature  measurements  shown  in  Fig. 
11  indicate  large  spatial  variations,  which  certainly 
exist  along  the  measurement  paths  as  well.  Numer¬ 
ous  researchers  have  attempted  to  make  corrections 
for  boundary  layer  effects.40’4143’44  However,  these 
techniques  are  limited  and  require  prior  knowledge  of 
the  general  shape  of  the  spatial  temperature  profile 
to  be  effective.  The  problem  of  solving  for  a  temper¬ 
ature  distribution  along  a  line  of  sight  by  probing 


Fig.  12.  FFT  power  spectrum  (0-150  Hz  shown)  of  the  temper¬ 
ature  versus  time  data  plotted  in  Fig.  10.  Fuel/air  equivalence 
ratio:  =  1.02;  position  C. 
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Fig.  13.  Temperatures  calculated  assuming  a  uniform  path  and 
using  various  area  ratio  combinations.  Fuel/air  equivalence  ratio: 
=  1.02;  position  C. 


multiple  absorption  features  may  be  described  with  a 
matrix  equation44: 


Si(Ti)  S1(T2)  S1(71re) 

“area!  " 

£2(^1)  S2(T2)  S2(Tn) 

PcJ-“2 

= 

area2 

_Sm(T1)  S2(T2)  Sm(Tn)_ 

_P„Lri_ 

_aream_ 

Two  area  measurements  allow  for  the  solution  of  two 
measurement  parameters,  T1  and  PiL1?  where  Lx  is 
known  in  this  case  (total  path  length): 

rw) 

_*S2(T1) 

Dividing  the  first  row  by  the  second  row  in  Eq.  (6) 
yields  the  familiar  simplification,  used  in  ratio  ther¬ 
mometry,  to  infer  an  average  temperature  with  two 
lines.  In  this  case,  the  area  ratio  is  equivalent  to  line 
strength  ratio,  as  expressed  in  Eq.  (4). 

The  scramjet  sensor  utilized  in  this  work  provides 
three  integrated  absorbance  measurements  and  may 
be  solved  to  yield  three  parameters.  Assuming  that 
we  want  to  solve  for  two  temperatures,  Tl  and  T2 
certain  parameters  must  be  prescribed  to  solve  for  a 
third  parameter  amongst  the  unknowns  Pi,Ll9p2, 
and  L2: 


pxLx  = 


areax 

area2 


(6) 


^i(^i)  *St(T2) 

Pi  Li 

P>L-2 

"areax" 

^2(^1)  $2(^2) 

_®3(Ti)  S3(T2)_ 

— 

area2 

area3 

For  example,  the  path  lengths  of  each  temperature 
bin,  L1  and  L2,  may  be  exactly  prescribed  and  the 
assumption  made  that  concentration  is  a  constant, 
Pi  =  Pz-  Therefore,  the  remaining  three  terms  to  be 
solved  are  7\,  T2,  and  p1  =  p2  =  p .  In  reality,  for  a 
combustion  product  like  water,  species  partial  pres¬ 
sure  is  also  a  function  of  temperature,  as  evidenced 
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Fig.  14.  Numerical  solution  assuming  two  temperature  bins 
(T1  and  T2)  of  size  90%  and  10%  of  the  total  path  length  and  a 
uniform  water  concentration.  Fuel/air  equivalence  ratio: 

=  1.02;  position  C. 


by  the  correlation  between  these  quantities  in  Figs. 
9-12.  This  situation  results  in  an  underdetermined 
system  with  four  unknowns  and  three  equations.  One 
may  use  thermodynamic  arguments  to  assume  an  a 
priori  temperature  functionality  for  the  partial  pres¬ 
sure,  p{T),  but  this  results  in  an  overdetermined  sys¬ 
tem  with  two  unknowns  (T1  and  T2)  and  three 
equations. 

As  an  initial  step,  it  is  useful  to  assume  a  uniform 
temperature  distribution  along  the  measurement 
beam  path  in  the  scramjet  and  infer  temperatures  us¬ 
ing  the  1392-nm/1343-nm  line  pair,  the  1392-nm/ 
1469-nm  line  pair,  as  well  as  the  1343-nm/1469-nm 
line  pair.  An  example  result  is  shown  in  Fig.  13.  The 
temperatures  that  are  calculated  with  the  1469-nm 
line  exhibit  higher  average  temperatures  and  larger 
scatter  than  the  1392-nm/1343-nm  line  pair.  The 
increased  scatter  is  to  be  expected  based  on  the  re¬ 
sults  shown  in  Fig.  13,  which  indicated  reduced  SNR 
for  the  direct-absorption  measurements  at  1469  nm. 
It  is  reasonable  to  expect  that  the  1392-nm  line 
yields  the  lowest  temperatures  since  this  low-E"  line 
is  particularly  sensitive  to  lower  temperatures.  A  re¬ 
view  of  Fig.  6  shows  that  the  line  strength  of  the 
1392-nm  feature  increases  greatly  at  lower  temper¬ 
atures.  Compared  with  the  higher-^"  lines,  the  pres¬ 
ence  of  locally  cold  regions  disproportionately  affects 
the  1392-nm  line. 

By  assuming  that  the  measurement  path  contains 
a  two-temperature  distribution  of  prescribed  lengths 
(L1  and  L2),  and  that  the  water  concentration  is  con¬ 
stant  throughout  (p1  =  p2  =  p),  it  is  possible  to  solve 
the  nonlinear  system  of  equations  in  Eq.  (7)  for 
7\,  T2,  and  p.  The  data  used  to  create  Fig.  13  is 
solved  for  these  quantities  using  a  Levenberg- 
Marquardt  nonlinear  routine,  as  shown  in  Fig.  14. 
The  simple  model  used  to  solve  these  equations  as¬ 
sumes  that  90%  of  the  path  length  is  at  one  temper- 


Fig.  15.  Temperature  versus  time  using  scanned-2/1  wavelength 
modulation  spectroscopy.  The  2 f  peak  ratios  of  water  features  at 
1392  and  1469  nm  are  used.  Fuel/air  equivalence  ratio: 

=  1.02;  1-kHz  bandwidth;  position  C. 
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Fig.  16.  SNR  comparison  for  direct  absorption  area  and  2 f  peak 
height  measurements  as  a  function  of  peak  absorbance  magnitude. 
Fuel/air  equivalence  ratio:  =  1.02;  1-kHz  bandwidth;  position  C. 


ature,  7\,  and  that  10%  of  the  path  is  at  T2.  This 
would  reflect,  for  example,  the  presence  of  a  hot 
boundary  layer  next  to  the  combustor  side  walls.  Nu¬ 
merical  fluid  simulations  of  the  scramjet  combustor 
predict  that  the  fuel  from  the  cavity-based  injectors 
tends  to  migrate  toward  the  side  walls  of  the  duct. 
This  would  lead  to  greater  heat  release  at  the  walls, 
resulting  in  higher  temperatures,  as  indicated  by  this 
idealized  interpretation  of  the  absorption  data.  In 
addition,  deceleration  of  the  flow  in  the  side-wall 
boundary  layer  would  cause  kinetic  energy  to  be  con¬ 
verted  into  increased  thermal  energy  (and  thus 
higher  temperatures).  By  probing  only  three  water 
features,  it  is  difficult  to  perform  a  realistic  analysis 
of  line-of-sight  distributions.  Although  the  solution 
used  for  Fig.  14  is  numerically  efficient,  given  the  few 
data  points  available,  an  accurate  result  clearly 
would  require  the  use  of  additional  lines.  These  lines 
should  be  selected  to  provide  good  SNR  as  well  as 
sensitivity  to  the  nonuniform  temperature  field.  Ad¬ 
dressing  these  issues  is  the  subject  of  future  work. 
Nonetheless,  the  current  data  set  provides  some  in¬ 
dication  of  higher  temperatures  near  the  walls,  as 
well  as  serves  to  illustrate  the  potential  of  multiwave¬ 
length  absorption  techniques. 

D.  2 f  Spectroscopy 

Although  the  large  absorbance  levels  in  this  rela¬ 
tively  long-path  combustor  make  direct-absorption 
scans  an  ideal  choice  for  obtaining  temperature  and 
concentration  measurements  for  the  strongest  wa¬ 
ter  absorption  features,  a  few  runs  were  analyzed 
using  scanned-2/*  wavelength  modulation  spectros¬ 
copy  to  verify  the  benefits  of  this  strategy  in  the 
case  of  weak  absorption,  as  demonstrated  by  the 
1469-nm  line.  Figure  15  shows  a  quarter-second  of 
time-resolved  2 f  temperature  measurements  using 
the  1392-nm/1469-nm  line  pair,  filtered  from  a 
2-kHz  scan  rate  to  a  1-kHz  measurement  band¬ 


width  that  is  analogous  to  the  direct-absorption  plot 
of  Fig.  9.  The  measurement  location  and  the  combus¬ 
tor  conditions  for  the  measurements  shown  in  Figs.  9 
and  15  are  identical,  and  the  average  temperature 
measured  in  Fig.  9  (direct  absorption)  is  used  to  cal¬ 
ibrate  the  2 f  measurement.  The  standard  deviation 
variation  in  the  average  temperature  in  Fig.  15  is 
21  K,  which  is  less  than  the  corresponding  variation 
of  40  K  with  the  direct-absorption  results  shown  in 
Fig.  9. 

The  true  strength  of  the  2f  method  lies  in  its  ability 
to  significantly  improve  SNR  for  weakly  absorbing 
features.3’24  As  detailed  in  an  earlier  work,14  2/*  peak 
height  ratios  are  simply  interpreted  as  being  propor¬ 
tional  to  line  strength  ratios  if  modulation  depths  are 
optimized  AND  if  absorption  levels  are  low  (less  than 
10%  absorption  typically).  Therefore,  the  long-path 
scramjet  combustor,  in  which  absorption  levels  are 
relatively  high  for  the  lines  utilized  in  this  study,  is 
not  an  ideal  test  bed  for  ratio  thermometry  with  2 f 
spectroscopy.  However,  by  comparing  the  SNR  of 
direct-absorption  measurements  versus  2 f  peak- 
height  measurements,  for  the  strong  lines  at 
1392  and  1343  nm,  as  well  as  for  the  weaker  line  at 
1469-nm,  it  is  apparent  that  the  2 f  technique  is  su¬ 
perior  for  weaker  features.  Figure  16  shows  a  com¬ 
parison  of  the  SNR  (spectrally  integrated  absorbance 
or  2 f  peak  heights  as  a  function  of  time)  of  the  indi¬ 
vidual  water  features  for  the  identical  combustion 
runs  shown  in  Figs.  9  (direct  absorption)  and  15  (2 f 
spectroscopy).  For  this  combustor  condition,  the  peak 
absorbance  of  the  1392  and  the  1343-nm  feature 
are  comparable  (—0.3).  The  SNR  of  the  spectrally 
integrated  absorption  strengths  (area  under  the 
Voigt  curve)  of  many  measurements  over  time  is  sim¬ 
ilar  to  the  SNR  of  2/*  peak  heights  over  time  (—22). 
For  the  weaker  feature  at  1469  nm,  the  direct- 
absorption  SNR  is  much  degraded,  whereas  the  2 f 
SNR  is  unchanged  (hence,  the  2 f  temperature  plot  in 
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Fig.  15  utilizes  the  1469-nm  line).  These  results  pro¬ 
vide  evidence  that  in  realistic  combustion  environ¬ 
ments,  2f  methods  improve  SNR  when  one  is  dealing 
with  weak  absorption.  The  high-frequency  phase- 
sensitive  detection  utilized  in  modulation  spectros¬ 
copy  helps  to  reject  noise.  Additionally,  and  perhaps 
more  important,  the  derivative  nature  of  2 f  spectros¬ 
copy  eliminates  the  need  to  perform  baseline  fits, 
which  is  a  large  source  of  uncertainty  when  one  is 
dealing  with  weak  transitions  or  broadened  spectral 
features.14’45 

The  scanned-2/*  method  used  here  requires  calibra¬ 
tion  and  is  a  relative  temperature  measurement, 
whereas  direct-absorption  scans  yield  absolute  tem¬ 
perature  as  well  as  species  concentration  informa¬ 
tion.  With  some  additional  work,  it  should  be  possible 
to  make  absolute  temperature  and  concentration 
measurements  with  scanned-2/*  spectroscopy.  If  ex¬ 
periments  are  conducted  to  validate  or  obtain  accu¬ 
rate  spectroscopic  parameters  such  as  broadening 
widths,  exponents  for  the  temperature  dependency  of 
the  broadening  widths,  pressure  shift  parameters, 
lower-state  energies,  and  line  strengths,  it  is  possible 
to  simulate  2 f  line  shapes  for  arbitrary  conditions, 
provided  that  the  exact  wavelength  modulation  depth 
and  intensity  modulation  depths  are  recorded.  Since 
2 f  signal  strengths  are  proportional  to  detected  laser 
intensity,  2 f  signals  must  be  normalized  by  the  dc 
signal  strength,  or  equivalently,  in  many  cases,  the  If 
signal.34  A  single  calibration  measurement,  prefera¬ 
bly  at  a  condition  similar  to  the  measurement  condi¬ 
tions,  is  all  that  is  necessary  to  establish  a  scaling 
factor  between  2f  measurements  and  2f  simulations. 
This  is  similar  to  what  is  done  for  large-modulation- 
depth  2 f  spectroscopy  of  broad  spectral  features,  as 
detailed  in  a  previous  article.44 

5.  Conclusions 

A  tunable-diode-laser-based  sensor  has  been  devel¬ 
oped  for  line-of-sight  temperature  and  concentra¬ 
tion  measurements  of  water  vapor  in  a  model 
scramjet  combustor.  Both  direct  absorption  scans  at 
4  kHz  and  2f  wavelength  modulation  scans  at  2  kHz 
were  used  to  obtain  line-shape  profiles  of  three  iso¬ 
lated  water  features  near  1392,  1343,  and  1469  nm. 
Integrated  direct-absorbance  area  ratios  were  used  to 
infer  temperature  and  water  vapor  concentration  at  a 
variety  of  engine  conditions  as  well  as  at  various 
vertical  path  locations  across  the  postcombustion  re¬ 
gion.  Optimized  optical  components  and  hardware 
have  yielded  excellent  signal  quality  in  the  probe 
beams  pitched  across  the  supersonic,  long-path- 
length  flow  field.  The  use  of  three  multiplexed  lasers 
in  the  current  sensor  has  allowed  for  a  preliminary 
discussion  of  line-of-sight  temperature  distributions. 
Results  from  scanned-2/*  measurements  indicate  im¬ 
proved  SNR  for  weakly  absorbing  transitions,  al¬ 
though  calibration  is  necessary  for  absolute 
temperature  and  concentration  measurements. 
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